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Abstract 

In this paper, one-parameter planar motion in generalised complex 
plane (or p—complex plane) Cp = {x + iy : x,y = p} which is 

defined as a system of generalised complex numbers is studied. Firstly, 
generalised Bobillier formula is obtained by using the geometric inter¬ 
pretation of generalised Euler-Savary formula in the p—complex plane. 
Moreover, it is shown that the Bobillier formula may be obtained by 
an alternative method without the use of Euler-Savary formula in the 
generalised complex plane. Thus, this formula generalises the complex, 
hyperbolic and Galilean cases. 
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1 Introduction 

A mechanism is a system of rigid elements (linkages) arranged and connected to 
transmit motion or is a combination of links which can transform a determined 
motion. In kinematics, a mechanism is a means of transmitting, controlling, or 
constraining relative movement. A planar mechanism is a mechanical system 
that is constrained so the trajectories of points in all the bodies of the system 
lie on planes parallel to a grand plane. The rotational axes of hinged joints that 
connect the bodies in the system are perpendicular to this ground plane, [T]. 

Kinematics concerned with the characteristics of movement without consid¬ 
ering the concepts of mass and force is sub-branch of mechanics and analyses 
the displacement of one point or a point system (object) with respect to time. 
Robot kinematics refers the analytical study of the motion of a robot manip¬ 
ulator. For a robot mechanism is very crucial for analysing the behaviour of 
industrial manipulators. 

Robot kinematics studies the relationship between the dimensions and con¬ 
nectivity of kinematic chains and the position, velocity and acceleration of each 
of the links in the robotic system, in order to plan and control movement and 
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to compute actuator forces and torques. The relationship between mass and 
inertia properties, motion, and the associated forces and torques is studied as 
part of robot dynamics, mia. 

The Euler-Savary formula which is one of the important kinematics formulas 
and has an important place in many fields such as mathematics, engineering 
and astronomy is found by Euler in 1765 and Savary in 1845, m- This formula 
giving the relationship between the curvatures of trajectory curves drawn by the 
points of the moving plane in the fixed plane is studied by [sjEjiiiiniiii] in two 
and three dimensional Euclidean spaces. Moreover, the Euler-Savary formula 
has too many applications in kinematics including in mechanical engineering of 
robotics in science of machines and mechanisms. 

In 1988, Fayet defined a formula giving the relation of the curvatures of 
second order of one-parameter planar motion and generalizing the Euler-Savary 
formula. Since this formula analytically solves the problem that the Bobillier’s 
construction solved graphically it is called the Bobillier formula, [T5] . 

In [13], it was demonstrated that the Bobillier formula can be obtained 
without the use of the Euler-Savary formula by Fayet. In addition, in m the 
Bobillier formula was shown by conventional procedures. 

Ersoy and Bayrak studied the Bobillier formula for the one-parameter planar 
motion in the complex plane, [5]. Also, they investigated the Bobillier formula 
in Lorentzian sense and saw that the same results can be achieved without the 
use of the Euler-Savary formula for the Bobillier formula. In doing so, they 
considered that the cases of polar curves to be timelike or spacelike, separately. 

Eri§ir et al. studied the one-parameter planar motion in the p—complex 
plane. So, they defined a one-parameter planar motion in the p—complex plane 
by transformation x' = [x — u) e*®'’ and obtained the relation between the 
velocities of any point in the moving p—complex plane as 

V“ = V^+V^ (1.1) 

where V“, and V)) are the absolute, the sliding and the relative veloc¬ 
ity vectors, respectively, [5D]. Moreover, in the generalised Euler-Savary 
formula giving the relation between the curvatures of the trajectory curves of 
one-parameter planar motion in the p—complex plane was given. 

In this study, we have obtained the generalised Bobillier formula by the 
geometric interpretation of the generalised Euler-Savary formula given in [21] 
for p—complex plane. In addition, we have expressed another presentation of 
the generalised Bobillier formula without using the Euler-Savary formula. The 
p—complex plane corresponds to complex, Galilean and hyperbolic plane for 
the special cases ofp= —1,0,-1-1, respectively. Since the Bobillier formulas in 
Galilean and hyperbolic plane have not existed according to our best knowl¬ 
edges, firstly we obtained them. Then we proved that the generalised Bobillier 
formula in p—complex plane satisfies the complex, Galilean and hyperbolic cases. 
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2 Preliminaries 


The generalised complex numbers or binary numbers are introduced as follows 

z = x + iy {x,y€R), i^ = iq + p {q,p€R). 

The double, dual and ordinary numbers are the particular members of two pa¬ 
rameter family of complex number systems. Moreover, the generalised complex 
number systems are isomorphic to the double, ordinary and dual complex num¬ 
bers when p + q'^/A is positive, negative and zero, respectively, [1]. 

Unless otherwise stated we assume that = p and q = 0 (p S R). This 
complex number system is denoted by 

Cp = {x + iy : x,y = p} . 

The set Cp is called the p—complex plane. For zi = (xi -I- iyi), Z 2 = (x 2 + 
* 2 / 2 ) S Cp the addition, subtraction and product are defined by 

Zi±Z 2 = (xi + iyi) ± {X 2 + m) =xi±X 2 + i{yi ± 2 / 2 )- 

and 

MP{zi,Z 2) = ixiX2 + P2/1P2) + iixiy2 + X2yi), 

respectively. The product definition yields the ordinary. Study and Clifford 
products as p is equal to —1, 0 and 1 

M~^{zi,Z 2 ) = ixiX 2 - yiy 2 )+ iixiy 2 +X 2 yi) for p = -l, 

M°{zi,Z 2 ) = ixiX 2 )+iixiy 2 +X 2 yi) for p = 0, 

M^{zi,Z 2 ) = {xiX 2 + yiy 2 ) + i{xiy 2 + X 2 yi) for p = l, 

respectively, [4]. For Zi = {xi,yi), 7.2 = ( 3 ^ 2 , P 2 ) G Cp, the scalar product is as 
follows 

MP(zi,^) = MP(^,Z2) = a;ia;2 -p2/i2/2 (2.1) 

and the magnitude of z = (x, y) £ Cp is given by the non-negative real number 

l|z||p = V|MP(z,z)| = -pp2| (2.2) 

where the over bar denotes the complex conjugation. The equalities 

|MP(zi,Z 2 )| = ||zi||p||z2||pCOsp0p, (2.3) 

||ziApZ2||p = ^P^||zi||p||z2|lpSinp6»p (2.4) 

are valid in p—complex plane where 9p is the p—rotation angle between the 
vectors Zi and Z 2 . 

In the Euclidean plane, a circle is defined in two ways. The first definition 
is the geometric locus of points which are at a fixed distance from a point. 
The second is the set of points from which a segment AB is seen at a constant 
directed angle 9p. Using these definitions, since a unit circle is the geometric 
locus of points z satisfying \\z\\p = 1, the unit circles in the p—complex planes 
are drawn as in Figure 1. 
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Let us interpret Figure 1. If we take p < 0, then the unit circle corresponds to 
unit ellipse of the form x^ + |p| = 1 and the complex number system Cp (p < 0) 
corresponds to elliptical complex number system. In this case, if p = —I, the 
unit circle in the p—complex plane corresponds to a standard unit circle defined 
by = I and the p—complex plane corresponds to the complex plane. 

If we consider p = 0, the above circle definitions yield different sets of points 
in Galilean geometry. According to the first definition, the equation ||z||g = x^ 
is hold and the unit circle is as x = ±1, (Figure 1) and for the first circle we 
cannot talk of curvature. According to the second definiton, the set of points 
is called a cycle. Figure 2 illustrates the cycle (inflection cycle, red lines) in 
Galilean plane and this cycle is the (Euclidean) parabola, [23]. Therefore, the 
space Co is the parabolic complex number system and the p—complex plane 
corresponds to the Galilean plane. Moreover, as can be seen from Figure 1 
the parabolic complex plane is divided into two area by the imaginary axis. 
Finally, when p > 0, unit circles in Cp (p > 0) are the hyperbolas of the form 
|x^ — py^l = 1 whose asymptotes are y = ±x/y^ (red dashed line in Figure 1). 
Thus, the space Cp (p > 0) corresponds to hyperbolic complex number system. 
Especially, if we take p = 1, the p—complex plane is referred to as a well-known 
hyperbolic plane and the asymptotes of the unit circles separate the hyperbolic 
planes into four regions, [lO] . 

p—trigonometric functions ( p—cosine (cosp), p—sine (sinp ) and p—tangent 
(tanp )) are, respectively, defined by 


cos p9p = 


sinp0p = 

tan pOp = 


cos (0p^/\p\) > 

1 , 

cosh {0p^/p) , 

Op, 

^ sinh (Op^) , 

ipOp 


cos p9p 


p < 0 

p = 0 {branchI) ; 
p>0 {branch I) 

, p < 0 

p = 0 {branch I) > 
p>0 {branch I) 


(2.5) 
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Additionally, the derivatives of p—trigonometric functions are as follows, 

(cos p6»p) =p sin p6>p, -^(sinp0p) = cosp6»p (2.6) 

dwp dup 

and the generalised Eulers formula is as follows 

= cosp9p + isinpOp. (2.7) 

The polar and exponential forms of any p—complex number z are 
z = rp[cosp9p + isinp0p) = rpC*®*’ 

where Vp = ||z||p and 9p are the p—magnitude and argument of z, respectively, 

m- 


3 The Generalised Euler-Savary Formula in p-Complex 
Plane 


Let T and T' be the moving and fixedp—complex planes, respectively. {O; ti, t 2 } 
and t' 2 } are the perpendicular coordinate system of these planes, re¬ 

spectively. Let A^p, Np and be the fixed points in the moving p—complex 
plane T and 7 p, 7 p and 7 p be centres of curvature of trajectory drawn by these 
hxed points in the fixed p—complex plane T'. 

The normals of trajectories drawn by these points pass from the instanta¬ 
neous rotation centre I and there is a p—rotation pole at each t moment. The 
geometric loci of the pole points in the fixed p—complex plane T' and the mov¬ 
ing p—complex plane T are called fixed pole curve {Q') and moving pole curve 
(Q), respectively, [20] . 

The pole curves rolls upon each other without sliding during the motion of 
p—complex planes T, T'. So, the pole curves (Q), {Q') are tangent to each other 
and have the same velocity at each t moment, [5D]. Thus, the real axis is the 
common tangent and the imaginary axis is the common normal for the pole 
curves. 

If 9p is the p—rotation angle of motion of the p—complex plane T with respect 
to T' at each t moment, then each point makes a rotation motion with 9p 
angular velocity at the instantaneous centre I. 

Let 


K = 


IN^ 

XJ. 

IlNlI 


K = 


XJ. -^p 


x^ = 


IN"' 

XX-^p 

|in! 


(3.1) 


be the unit vectors in the direction of the pole rays IN^, IN^ and IN^, respec- 
tively, (Figure 2). So the exponential forms of these vectors are given by 


X 


1 

p 


le 


iOt 


= ie^ 


iSt 


X^ = ie 




(3.2) 


where = p and 9^, 0^, 9p, are p—rotation angles of motion in the p—complex 
plane. 
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The p—distances of the points and 7 ^ from the origin I are pk and 
{k = 1, 2, 3), respectively then from the equations (2.2) and (2.3) we can write 


Similarly 


and 


MP 

(iNj.Xj) 

= Pi, 

MP ( 

i7p\x;,) 

MP 

(iNj.xj) 

= P2, 

MP 


MP 

(INJ.XJ) 

= P3, 

MP ( 



P'l- 


= P 2 


= Ps- 


(3.3) 


An inflection point may be defined to be a point whose trajectory momen¬ 
tarily has an infinite radius of curvature. Such points also have zero acceleration 
normal to their trajectory, [ZllIHllll]. Let the inflection points be N^*, N^* and 
Np*, (Figure 2). The locus of such points is a circle in the moving p—complex 
plane T called as an inflection circle. So, from the equation (2.2) and (2.3) we 
can write 


|m^(in;*.Xi)| =p*. |mp (iNf ,X2)| =p*, |m^ (inJ*.X 3)| =p; (3.4) 

for the abscissae of the inflection points. 
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Figure 2. The inflection circles of motions in the generalised complex plane. 


There is a relationship between h which is diameter of the inflection circle 
and p* (the p—distance from the instantaneous centre I to inflection point ). 
Namely; the sliding velocity vector of the inflection point N^* is perpendicular 
to the vector which connects the centre to this point. Moreover, the sliding 
velocity vector is perpendicular to the angular velocity vector. So, we can write 

= Wpple^^l (3.5) 

where Wp = WpZp is angular velocity vector of motion in the p—complex plane 
and Zp is the unit vector in the direction of the angular velocity vector. 

On the other hand the points move around an trajectory whose instantaneous 
centre is I during the motion T/T' in the p—complex plane. So, the amount of 
displacement of the point I is equal to the product of the diameter h with the 
amount of angular displacement. Thus, we can write 

Vj^ = hwp. (3.6) 

Considering the equations (3.5) and (3.6), the relation between p* and h is 

h = p*e*®p. (3.7) 

During one-parameter planar motion T jT' in the p—complex plane, the point 
Np in the moving p—complex plane T draws a trajectory which has curvature 
centre 7 ^ in the fixed p—complex plane T'. Conversely, the point 7 ^ in T' draws 
a trajectory with curvature centre Np in T. This relation between the points 
Np and 7 ^ is given by 


1 

Pi 



le 



(3.8) 


where R and R' are the radii of curvature of the pole curves {Q) and (Q')j [21] • 
Considering the equation (3.6), (3.7) and (3.8) the following theorem can be 
given. 

Theorem: Let h be the diameter of the inflection circle in p—complex plane 
and Wp be the angular velocity of motion in the p—complex plane. The first 
form and second form of the generalised Euler-Savary formula are 





1 

R 


and 


1 

h 


Vi 


respectively. 

This theorem can be interpreted as the following special cases; 


(3.9) 


(3.10) 
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Case 1. We consider p = —1. In that case we reduce that the generalised 


VPi PiJ 


le = -k — 


R' 


Euler-Savary formula given in the equation (3.8) to 

where i is the complex imaginary unit satisfying = —1. So, the generalised 
Euler-Savary formula coincides with the Euler-Savary formula in the complex 
plane, m- Moreover, from the equations (3.9) and (3.10) it is seen that the 
first form and second form are ^ ^ and ^ as given in |16) . 

Case 2. Let p = 0. Erom the equation (3.8), we obtain = 

e (-^ — ■^) , where e is the dual imaginary unit satisfying = 0 which cor¬ 
responds to the Euler-Savary formula in Gallilean plane, [9]. Moreover, con¬ 
sidering the equations (3.9) and (3.10) the first form and second form are 
I = and ^ = y- In addition, we calculated the first form and 

second form of Shear motion in Galilean plane on aspect of [3] and we saw that 
the above equations are the same with calculated equations by us. 

Case 3. If we take p = 1, then we obtain that the Euler-Savary formula is as 




~ W ~ 'k j is the hyperbolic imaginary unit satisfying 

= 1, j ^ ±1. This formula corresponds to the Euler-Savary formula under 
the planar motion in the hyperbolic plane, [8]. Also, from the equations (3.9) 
and (3.10), the first form and second form are equal to and ^ = y- 

Similar to special case 2, considering the reference [5], we obtained that the first 
form and second form for the motion in the hyperbolic plane and proved that 
they are coincident with the above equations. 


4 Generalised Bobillier Formula obtained by Gen¬ 
eralised Euler-Savary Formula 

We consider that the inflection points N^* and N^* are on the direction 
of (J, Xp), (/, Xp) and (/, Xp), respectively. So, the images of these inflection 
points are Qp, Qp and Qp where IQ* = ^X* , (1 < fc < 3), (Figure 3,4,5). 

Then, from the equation (2.2) and (2.3) we have 








M” 




(4.1) 


So, from the equation (3.7) the following equations hold; 




^ = iK 


= ^Xpc 

IQ^pe-< = = iX^ 

= |x^e 


1 X^e"*®p = iXt 


h^p- 


(4.2) 


The last three equations show that 




iQp.xi 


'’p = M’ 




> = M’ 


(iQp.x3)| 


el ^ ]_ 
h 


(4.3) 


So, the set of the points Qp is a straight line D parallel to axis x. Thus, the 
line D is an image of the inflection circle of the p—complex plane by means of 
an inversion at the rotation centre /, (Figure 3,4,5). 

Since the points Qp, Qp, Qp are linear, the vectors (IQi - IQ 2 ) and (IQ^ - IQ^) 
are linearly dependent. So, the p—cross product of these vectors Is 

(IQi - IQl) Ap (IQ2 - IQ3) = 0. 


Thus, we can write 

(lQ>pIQ2) + (iQ^pIQi) + (iQ^pIQ^) = 0. 

If we consider that P 1 P 2 P 3 ^ 0, since IQp = IQp = and 

IQp = we can easily write 


Pi 


+ P2 


„3 

P A„ie 


+ P3 


Considering the last equation and the equation (2.4) the following theorem can 
be given. 

Theorem: In one parameter motion T /T' in the p—complex plane the 
relation between the centres of curvatures Is 


p\ slnp 0 . 


|23 

p 


■ P *2 sinp^^^ + P 3 sinpef = 0 


(4.4) 


where ^ - ]- and are the p—rotation angles between XL and X™ 

for the permutations of the indices fc, Z, m = I, 2,3 ; 2,3,1; 3,1, 2. 

The formula (4.4) Is called of the generalised Bobilher formula which is 
generalisation the generalised Euler-Savary formula for all p G R. The special 
cases of the generalised Bobillier formula with respect to the sign of real number 
p are as follows; 

Case 1. If we take p < 0, the p—complex number system Cp corresponds 
to the elliptical complex number system. From the equation ( 6 ), since sinp0p = 

sin (J)p^/\j^\j the generalised Bobilher formula is equal to 

pI slnp (of^Ap\^ + p *2 sinp + pi sinp = 0 - 

Especially, if we take p = — 1, this formula becomes 

p* sln0^^ + P 2 sln0^^ + pI sln0^^ = 0. 


Under the circumstances, the generalised Bobillier formula is reduced to the 
Bobillier formula at the complex plane. Figure 3 represents the Illustration of 
the Bobillier construction in the complex plane, |16] . 
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Case 2. If we consider p = 0, the p—complex number system Cp is equal 
to the parabolic complex number system. By considering the equation (2.5) if 
we substitute sinp0p = 0p into the equation (4.4), we get 

+ +P3^p'=0. 

This formula is the Bobillier formula for Shear motion in Galilean plane. The 
inflection cycle for Shear motion is illustrated as 



Case 3. When p > 0, the p—complex number system Cp is referred to as 
the hyperbolic complex number system. In addition, considering the equation 
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(2.5), we have sinpOp = sinh (6>pY/p). Thus, we obtain that the generalised 
Bobillier formula is reduced to 


p\ sinh (Ol^y/p) + p *2 sinh {Of^) + pi sinh = 0. 

Especially, if we take p = 1, we obtain that the Bobillier formula as 

p\ sinh + P 2 sinh 9^^ + pi sinh 9^^ = 0. 

Thus, this formula is equal to the Bobillier formula at the hyperbolic plane. The 
figure of the (hyperbolic) inflection circle for p = 1 is as follows 



5 An Alternative Way for the Generalised Bo¬ 
billier Formula 

In this section, we will obtain the generalised Bobillier formula in the p—complex 
plane without the use of the generalised Euler-Savary formula. So, the gener¬ 
alised Bobillier formula will be obtained considering the velocities and acceler¬ 
ations of planar motion in the p—complex plane. 

Firstly, we calculate the trajectory velocities and accelerations of the points 
in the moving p—complex plane T. Let V“ (./Vp) and Jp (iVp) be absolute 
velocity and acceleration vectors of the point iVp, respectively. Moreover, if Wp 

dO^ 

is the angular velocity of motion T jT' in the p—complex plane, Wp = where 
9p is p—rotation angle. So, considering the equation (3.5) the sliding velocity 
vector of the point Vp is 

(Vpi) = wpp.e^^p (5.1) 
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In addition, from the equation (1.1) there is the relation between the velocity 
vectors 

KiK)=%iIp)+^iiNp) ( 5 - 2 ) 

where V“, and V)) are the absolute, sliding and relative velocity vectors, 
respectively. If we consider the equations (5.1) and (5.2) it is obtained that 

V;{N^)=V;{Ip) + Wj,p,e<. (5.3) 

Differentiating the equation (5.3) with respect to t, we get 

j; (N^) = j; (Ip) + (5.4) 

where the first term is the path wise tangential acceleration component, the 

second term is the centripetal component and the third term can be shown to 
be a pure imaginary component. Considering this analysis, the absolute velocity 
vector of the inflection point is linearly dependent with the absolute acceleration 
vector of the inflection point since the normal component of acceleration is zero. 
So we can write 

v; (ivi*) Apj; (ivi*) = 0. (5.5) 

Here considering V)) = 0 for the inflection point Np* and using the equations 
(5.3) and (5.4) we can easily find 


Wppl ( e* P ApJp (Ip) ) + plwpWp ( e* P Ape' P ) + pI'^w^ ( e* P Apie' P J = 0 


If we consider the equation (2.4) and make necessary arrangements we obtain 
that 


WpPl (e<ApJp(/p)) + pfw. 
and finally it is said that 


pZp — 0 


PlZp — 


e^‘^Mpj;(/p) 

w'i 


(5.6) 


With similar process for the inflection points N'i* and iV^* the following equa¬ 


tions hold; 


P2Zp — 


e^'^pApj; (Ip) 


wt, 


(5.7) 


and 


P^Zp — 


e-^pApj;(/p) 


(5.8) 


On the other hand the linear connection between X^, and X^ may be written 
as follows 


Aie*^p -f Aae*''; -f Age*''? = 0 




iSt 
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where AijAsjAa G K. From the last equation by successive p—cross product 
with and we get 

Ai = sinp0p^, A2 = sinp0p^, A3 = sinp0p^ 

where 9^^ =9^ — 9^, 9^ = 9^ — 9^ and 9^ = 9^ — 9p. So, we obtain that 

sinp6lp^e*^*> + sinp6lp^e*®p + sinp0p^e*^p = 0 


The p—cross product of the last equation with is 


sinpe 


PApJ {Ip 


+ sin 2 


P ApJ (Ip) 


+ sin 2 


P ApJ (Ip) 


= 0 . 


Finally, if the equations (5.6), (5.7) and (5.8) are substituted into the last equa¬ 
tion, the following formula is hold; 

pI smp9p^ + P 2 sinp0p^ -b pg sinp0p^ = 0. (5.9) 

This is the formula given in (4.4). So, this direct way gives us the generalised 
Bobillier formula without using the generalised Euler-Savary formula in the 
p—complex plane. 


6 Conclusion 

The angle between the tangent of pole curve at the instantaneous pole centre of 
coupler with respect to the base of a four bar linkage and of the cranks is equal 
to the angle between the other crank and the collineation axis. This expression 
is can be verified by the Bobilliers construction by graphically. This has been 
a major interest to physicists, mathematicians and engineers. Also this prob¬ 
lem can be solved by an analytical method called the Bobillier formula which 
is more practical to use. Various geometric and analytical methods have been 
developed for the Bobillier formula of Euclidean and Lorentzian planar motion. 
In literature, the Bobillier formula in Galilean and hyperbolic planes has not 
been obtained yet. Then the following questions can be asked. ’’Can be this 
formula obtained in Galilean and hyperbolic planes?” and ”Is it possible to give 
a generalised formula for the Bobilliers construction in a p—complex plane in¬ 
cluding all planes?” Thus, we find out the Bobillier formula for p = 0 (Galilean) 
and p = 1 (hyperbolic). Moreover, since the p—complex plane correspond to 
complex, Galilean and hyperbolic planes for the special cases ofp (—1, 0, 4-1), re¬ 
spectively, we obtain that the generalised Bobillier formula in p—complex plane. 
We give the special cases of the generalised Bobillier formula for the cases of 
p. So, we think that this study would be important to the sciences of mathe¬ 
matics, engineering and astronomy. There is no doubt that this study will add 
innovation to science. 
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